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Abstract 

The density of complex eigenvalues of random asymmetric NxN matri- 
ces is found in the large- ./V limit. The matrices are of the form Hq+A where 
A is a matrix of -/V 2 independent, identically distributed random vari- 
ables with zero mean and variance N^ 1 v 2 . The limiting density p(z,z*) 
is bounded. The area of the support of p(z,z*) cannot be less than nv 2 . 
In the case of H commuting with its conjugate, p(z,z*) is expressed in 
terms of the eigenvalue distribution of the non-perturbed part H Q . 
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Random Hermitean and real symmetric matrices have been extensively stud- 
ied since the 50's, the time when Wigner introduced them into theoretical 
physics. A lot of results concerning these matrices and respective techniques 
are known now. In contrary to this, random complex and real asymmetric ma- 
trices are much less studied. Although they have already proved to be useful. 
We mention here only two examples (but see a discussion in jl]]). These are: i) 
quantum chaotic scattering and decaying processes, where complex eigenvalues 
of random non-Hermitean matrices are used to analyse statistical properties of 
resonances jl], 0, B|, and ii) neural network dynamics where synaptic matrices 
are in general asymmetric and the distribution of their eigenvalues is important 
for the understanding of network dynamics j|, ^] . 

In this letter we consider random real asymmetric matrices of the form H = 
Hq + A. A = [a>jk]jk=i is a matrix of iV 2 independent, identically distributed 
random variables such that 

(a jk ) = 0, (a jk ai m ) = N~ 1 v 2 5 j i5 km . (1) 

The angle brackets (. . .) denote average over the random variables aj k . For 
simplicity we assume that aj k are Gaussian but our results remain valid for a 
wider class of distributions. We treat A as a perturbation and Hq as the non- 
perturbed part and our aim is to determine the large- N limit of the averaged 
density of complex eigenvalues of Hq + A. 

If A and Ho are symmetric (or Hermitean) and A obeys the GOE (GUE) 
statistics, then Hq + A is known as the deformed GOE (GUE) ||. In this 
ensemble eigenvalues are real and, hence, their density is completely determined 
by the imaginary part of the Green's function G(E + iO) = (iV~ 1 tr [E + iO — 
H)^ 1 ). It appears that in the large- N limit the Green's function of the deformed 
GOE (GUE) is related to that of Hq by the so-called Pastur's equation |7j: 

G(z) =G (z-v 2 G(z)). (2) 

Although this equation cannot be solved explicitly (except of a few cases) it 
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provides useful information about the density of eigenvalues. For instance, 
one can prove that the density of eigenvalues in the deformed GOE (GUE) 
is bounded and generically decays as the square root in the vicinity of the 
spectrum boundaries ||. 

Eigenvalues of asymmetric matrices are complex and their average density 
p(z,z*) is determined by the electrostatic potential 

$(k,z,z*) = - N- 1 {log det[(z I - H)*(z I - H) + k 2 I}) 

by means of Poisson's equation p(z,z*) = [0> 0' ^ * s 

identity matrix. Positive infinitesimal k is introduced in order to regularize 
the potential. Provided k = 0, $ as a function of complex z has a singularity 
whenever z equals one of the eigenvalues of H. 

Anticipating an important role of positive semi-definite matrices 7i = (zl — 
H)*(zl — H) in studying complex eigenvalues of H we introduce the following 
Green's function 

R{k) = (N-hr (H + k 2 /) -1 ) (3) 

corresponding to Ti. R(k) as a function of k is analytic in the right half 
of the complex plane and obviously determines the density of eigenvalues of 
TL. We show that in the large-iV limit the Green's functions of TC and 7io = 
(zl — Hq)*(zI — Hq) are related by the equation ([!(]) which can be thought 
as generalization of Pastur's result to the case of positive semi-definite random 
matrices. In passing we find derivatives of the electrostatic potential. This al- 
lows us to derive an expression for the average density of complex eigenvalues of 
H, p(z, z*), and for the domain of their distribution. The respective expressions 
(|l^) - (|l4] ) are given in terms of Hq. Actually they set up the only restriction 



to Hq: quantities entering ( |i2| ) - (14) must be well defined in the large- N limit. 
We do not specify Hq further. It can be real or complex and either deterministic 
or random. In the latter case it is assumed that Hq is statistically independent 
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of A and it is understood that the average over realizations of Hq has been 
taken. 

At this point it is worth mentioning that in the specific case of Hq commuting 
with its conjugate Hq* (i. e. Hq can be symmetric, skew-symmetric, Hermitean, 
skew-Hermitean, etc) p(z, z*) can be explicitly expressed in terms of the density 
of eigenvalues of Hq (see (|l5|)-(|T7|)). This should be compared with the case of 
deformed GOE (GUE) where only the relation @ between Green's functions 
is known. 

Our last remark concerns matrices studied in ||. They are of the form 
iVV T + B, where V is an N x M matrix of NM independent Gaussian variables 
and B obeys the GOE statistics. These random matrices differ from those 
considered here in that aspect that B is symmetric while A is asymmetric. The 
eigenvalue distribution of VV T is known [1C] and it seems interesting to recover 



the results of [3| , which were obtained by means of the replica trick and 
supersymmetry calculations ||, in the framework of our approach. But this 
problem goes beyond the aim of the present letter. 

Introducing the notation Q{k) for the inverse to H + k 2 I we rewrite the 
following obvious matrix identity / = (G(k)(TL + k 2 /) -1 ) as 

k 2 (G(k)) =I-(zI- H )*((zI - H)Q(k)) + (A*(zl - H)G(k)) . (4) 

zl — Hq is statistically independent of G{n) but A which enters the (zl — H) 
term in the r.h.s. of (Q) is not. In order to decouple (AQ(k)) and (A*AQ(k)) 
we first notice that each of the entries Q pq of the matrix Q(k) is a function of 
the Gaussian variable ai m . Therefore 

(ai m Gp q ) = (af m )(dg pq /dai m ) = N~ l v 2 (dG pq /da lm ). (5) 

This is the only place where Gaussian distribution of a\ m is used. In the 
non-Gaussian case it can be shown that (||) holds up to the 1/iV 2 order if 
aim = N~ l l 2 a.i m and the random variables a\ m possess several first moments. 



5 



Straightforward application of (^) and the following rule for differentiating ma- 
trix elements of G(n) with respect to those of A 



dQ pq 
daim 

gives 



[G(zi - H)*] pl g mq + g pm [(zi - H)g] lg . (6) 



((zl - H)Q(k)) = (zl - H )(G(k)) - v 2 ((zl - H)Q{k) ■ N- L trg{K)) + 0(1/N). 

One can readily check (JgJ) making use of dQ pq / dTLkm = —GpkG mq and of the 
chain rule. The normalized trace of G(k) is a self- averaging extensive quantity 
Q. That is it becomes non-random in the large- N limit: iV -1 tr = R(k) + 
0(1/ N), where R(k) = (N'hr Q{k)). Therefore we conclude that 

((zi - H)G(k)) = (zl - H ){G(k))[1 + v 2 R(k)}- 1 + 0(1/N). (7) 

Similar reasoning leads to 

{A*(zl - H)G(k)} = -v 2 k 2 R(k){G(k)} + 0(1/N). (8) 

Collecting (|) and (0)-@ we find that in the leading order 

1 + v 2 R(k) 

{G{K)} = (zl - H )*(zl - H ) + k 2 [1 + v 2 R(k)] 2 I ' ( ' 

Introducing the notation Go(k) for the inverse to TLq + n 2 I one can write (||) 
in the form (£/(«)) = (1 + v 2 R(k))Gq(k[1 + v 2 R(k)\). R(k) is to be determined 
from the self-consistency equation 

R( K ) = [1 + v 2 R(k)]R (k[1 + v 2 R(k)]), (10) 

where R (k) = N-hvGo^)- 

Since -d^/dz* = (N~Ht(zI - H)G(k)) one can use (@) and © to calculate 
p(z, z*). Indeed, 

_ d$faz,z*) = N - hr ( zI _ Ho) g o{K[1 + V 2 R( K )]) + (1/N) 
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Simple analysis of ( |10D shows that in the leading order q z * — 

by 



K=0 



8$(k, z, z*) 



dz* 



K=0 



N-hr(zI-Ho)g ^(z,z*)), 



m given 



(11) 



where j(z, z*) = lim K ^ + K [l + v 2 R(k)] is the solution of Ro(sf) = v 2 if z lies 
inside the domain D determined by the inequality 



Ro(0) = N-hi[(zI - H Q )*(zI - H)}- 1 > v' 



(12) 



and "/(z.z*) = otherwise. Since in the latter case M% does not depend 

,v ' ' dz K =o 

on z we conclude immediately that p(z, z*) = outside D. On the other hand, 



differentiating (11) with respect to z one finds that inside D 

p(z,z*) = (ttv 2 )- 1 -*- x I{z,2?) 
where I(z,z*) is the large- N limit of 

N-hv (zi - H )g (j(z, z*))(zi - H yg ( 7 (z, z*)) 
n-\t{zI- H )g 2 ( 1 (z,z*)) 2 \N-hxg 2 ( 7 (z,z*)) 



(13) 



-1 



(14) 



For any two matrices P and Q |trPQ*| 2 < tr PP*trQQ*. Therefore I(z, z*) > 
and p(z,z*) is bounded by (ttv 2 ) -1 . This fact which is interesting in its own 
has also an important consequence: the area of D, the support of p(z, z*), is 
not less than the area of a disk with radius v. 

In order to illustrate the formulas derived we consider a few examples. If 
Hq = 0, then the equation i?o(0) = v~ 2 which determines the boundary of D 
takes the form \z\ 2 = v 2 and I(z,z*) obviously vanishes. Thus, we recover the 



circular distribution [11, |T^, [L3|]: p(z, z*) equals ( 



TTV 



2\-l ; 



inside the disk \z < v 



and zero outside. 

In our next example Hq is the Jordan block J = [/io(5j — i t fc]j v fc=1 , ho > 0. 
J has only one eigenvalue z = which is defective and highly sensitive to 
perturbations. On replacing zero in the left lower corner of J by small positive 
e, one gets iV distinct eigenvalues ^(e/ho) 1 ^ exp (2-Kik/N). For fixed N the 



perturbed eigenvalues approach zero as the parameter e/ho vanishes but the 
rate of convergence is extremely slow if N is large. For instance if N = 50 one 
needs e/ho oc 10 -50 in order to confine the eigenvalues into the disk |z|//io < 0.1. 
Therefore, if not exponentially small, perturbation splits zero eigenvalue of J 
into the circle \z\ = ho. As can be seen from (12) this phenomenon manifests 
itself in the large-iV limit. Indeed, in the case of Hq = J (|l2| ) reduces to 



\ z \ 2 ~ ^o 2 | < v<1 ■ Therefore if v < ho the eigenvalues of J + A are distributed 

2 2 

< 1 + -^-ty which degenerates into a circle as 
h 



in the annulus 1 







z 



v vanishes. When v > ho the eigenvalues are distributed in the disk 



z 

h~o~ 



< 



1 + 7^-7. In FigJT] we present results of numerical diagonalization of random 
h ^ 

matrices J + A. As can be seen, the correspondence between the numerical 
results and our analytical predictions (for N — > oo) is quite good. 
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Figure 1: Distribution of numerically computed eigenvalues of the random 
matrices J + A in the complex plane z/ho- In each of the plots N = 50 and the 
number of samples is 40. (a) v 2 jh\ = 1/2, (b) v 2 jh\ = 1. The full lines show 
the boundary of the support of p(z, z*) in the large- N limit. 



If Hq commutes with its conjugate Hq* , our formulas (p|)-(||) become 
simpler. Let us assume for certainty that the eigenvalues of Hq are real. Then 
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the boundary of D is determined by 

n(X)dX _ 1 



\z — A| 2 v 2 



(15) 



where n(X) is the density of eigenvalues of Hq. p(z,z*) is given by the same 
expression ( Jl3| ) as before but now I(z,z*) is 

\z-X\ 2 n(X)d\ 



I 



\z- A[ 2 + 7 2 (z,z*)] 2 

(z - X)n(X)dX 



\z- A| 2 + 7 2 (z,z*)] 2 
and 7(2, z*) has to be found from 



2 r 



n(X)dX 



1 -1 



\z - X\ 2 +7 2 (z,z*)] 2 



n(X)dX 1 

A|2 + 7 2 ~ y 2 



(16) 



(17) 
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